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' Abstract 

o 

D ' We present in this article a detailed quantitative discussion of the measurement of the leptonic 

, mixing angle 613 through currently scheduled reactor neutrino oscillation experiments. We thus focus 

on Double Chooz (Phase I & II), Daya Bay (Phase I & II) and RENO experiments. We perform a 
unified analysis, including systematics, backgrounds and accurate experimental setup in each case. Each 
identified systematic error and background impact has been assessed on experimental setups following 
^ ' published data when available and extrapolating from Double Chooz acquired knowledge otherwise. 

D I After reviewing the experiments, we present a new analysis of their sensitivities to sin^(2Si3) and study 

^ ^1 the impact of the different systematics based on the pulls approach. Through this generic statistical 

^ ' analysis we discuss the advantages and drawbacks of each experimental setup. 

C<\ '• 1 Experimental context 

> 
00 

' Over the last years the phenomenon of neutrino flavor conversion induced by nonzero neutrino masses 
T^lj- . has been demonstrated by experiments with solar [H O [3] , atmospheric [1] , reactor O E] and accelerator 
' neutrinos [71 \8\. Neutrino oscillation, that can be described by the Pontecorvo-Maki-Nakagawa-Sakata 
. (PMNS) mixing matrix [llj . is the current best mechanism to explain the data. Considering only the three 
' known families, the neutrino mixing matrix is parameterized by the three mixing angles {0i2, O23, ^13) and a 
. possible S CP violation phase. The angle 9i2 has been measured to be large (sin^(26'i2) — 0.8), the angle O23 
' has been measured to be close to maximum (sin^(2023) ^ 0.9); but the last angle 613 has only been upper 
bounded sin^(26'i3) < 0.15 at 90 % C.L., by the CHOOZ reactor experiment [9l[l2]. These achievements 
have now shifted the field of neutrino oscillation physics into a new era of precision measurements. Next 
generation experiments are underway all around the world to further pin down the values of the oscillation 
' parameters of both solar and atmosheric driven oscillations. Currently the most important task, for the 
experimentalists, is the determination of the last oscillation through the measurement of the unknown 
mixing angle ^13. An improved sensitivity on ^13 is not only important for the understanding of neutrino 
oscillations, but also to open up the possibility of observing CP-violation in the lepton sector if the ^13 driven 
oscillation is discovered by the forthcoming experiments. 

In order to improve the CHOOZ constraint on ^13 at least two identical unsegmented liquid scintillator 
neutrino detectors close to a nuclear power plant (NPP) are required. The near detector(s) located a few 
hundred meters away from the reactor cores monitor the unoscillated 17^ flux. The far detector(s) is(are) 
located at a distance between 1 and 2 km, searching for a departure from the behavior induced by 

oscillations. Experimental errors are being partially cancelled when using identical detectors. The goal is 
to achieve an overall effective systematic error of less than 1 % [T^ [31] . 

Three experiments have received partial or full approval to perform such a measurement in a near future: 
Double Chooz in France [21], Daya Bay in China [25] and RENO in Korea In addition, a project is 
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under study at the Angra power plant in Brazil to further improve the sensitivity of the measurement on a 
longer time scale [H]. Furthermore the Japanese KASKA collaboration is promoting the reactor neutrino 
oscillation physics [28] . 



2 Neutrino oscillation at reactor and ^13 



Fission reactors are prodigious sources of electron antineutrinos which have a continuous energy spectrum 
up to about 10 MeV. For > i?thr = 1-806 MeV they can be detected though the Ve+P e^+n reaction 
using the delayed coincidence technique, where an electron antineutrino interacts with a free proton in a 
tank containing a target volume filled with Gd loaded liquid scintillator. The positron and the resultant 
annihilation gamma-rays arc detected as a prompt signal while the neutron slows down and then thcrmalizes 
in the liquid scintillator before being captured by a hydrogen or gadolinium nucleus. The excited nucleus 
then emits gamma rays which are detected as the delayed signal. Electron antineutrinos energy is derived 
from the measured visible energy from positron energy loss and annihilation, 

E^i, = E,+ +me-E,- S*'"- + 2 me . (1) 



The Ve spectrum is calculated from measurements of the beta decay spectra of ^^^U, ^"^^Pu and ^^^Pu [T51 
after fissioning by thermal neutrons and theoretical ^ssy spectrum, since no data are available for ^"^^Ulll. 
As a nuclear reactor operates, the fission element proportions evolve in time (the so-called burn- up). Since 
wc are interested here on long term interpretation of the data for oscillation search, we will use an average 
fuel composition for a reactor cycle corresponding to 

235u (55.6 %), 23^Pu (32.6 %), '^^U (7.1 %) and ^^ipu (4.7 %) . (2) 

The mean energy release per fission, (Ef), is then 205 MeV and the energy weighted cross section for 
Vep — > e~^n amounts to (cr/) = 6 10"^'^ cm^ per fission. Let us introduce a new luminosity unit, called the 
r.n.u. (for reactor neutrino imit) and defined as 1 r.n.u. = 0.197 10^" MeV. With this imit, an experiment 
taking data for T years with a total NPP (nuclear power plant) thermal power of P GW and with N 10^° free 
protons inside the target has a luminosity C = T P N r.n.u.. The event number, N{L), at a distance L from 
the source, assuming no -oscillation, can be quickly assessed with 



For the full antineutrino reactor energy spectra simulation, we follow Vogel's analytical parameteriza- 
tion |16| . based on second order E^, polynomials. Higher order parameterizations [18| give very comparable 
results and do not require a specific attention for our aim in this article. The antineutrino event rates 
per energy range is then computed according to the mean reactor core composition ([2]) and experimental 
site specifications (reactor and detector locations, average efficiencies and running time as described in 
section [B]) . 

Reactor neutrino experiments measure the survival probability P^e which does not depend on the Dirac 
5cp phase. In addition, the oscillation of MeV's reactor neutrinos studied over a distance of a few kilome- 
ters is not affected by the modification of the coherent forward scattering from matter electrons [311 [33] . 
Expanding the full three flavors oscillation probability as a function of {Am2i/ Am'^^)'^ ~ 1/30^ ratio 
and sin^ (26*13), Pee measurements from reactor experiments on the kilometer scale may be described by the 
simplest two flavors oscillation formula: 

PeeC^l-sin2(2Msin2(^^) • (4) 

as long as sin^(26'i3) > 10~^. We assume that in cq. (d]), Arn^^ is measured by other experiments (MINOS [S], 
K2K [7] and super-K [J]). With a determination of Am^^ better than 10 %, the impact on sin^(20i3) 

1238^ fissions only with fast neutrons. Theoretical predictions arc computed by summing all known beta decay processes 
contributing. 
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determination can be neglected [31]. All results within this study are, unless otherwise stated, computed 
for a representative value of Ato|]^ [T^ : 

Am^i = 2.5t|^:25 10"^ eV^ at 68 % C.L.. (5) 



3 Generic analysis of 6*13 sensitivity 

The calculation of event rates is a convolution of the flux spectrum, the cross section, the oscillation 
probability, the detector efficiency with the energy response function. 

The detector energy response simulation, as well as its correction through detector callibrations are specific 
to each experiment. The details of the corrections are thus beyond the scope of this article. We therefore 
assume in the following that the reconstructed energy is identical to the true deposited energy. We will 
implement, anyway, a simple energy scale systematic uncertainty (section l6.ip . 

We based our event rates computations on an extended version of the numerical code developed for Dou- 
ble Chooz [311 [HI ■ These computations take into account the characteristics of each experimental setup as 
the number of reactors, detectors, locations, overburdens, efficiencies, operating time, and so on which will 
be described in section [6l 

The resulting event rates, then, form the basis for a analysis, where systematic uncertainties are properly 
included. Since event rates in the disappearance channel of reactor experiments are quite large, we can use 
a Gaussian x^, which has the advantage of allowing a natural inclusion of systematic errors through the 
so-called pull- approach" [H] : 



Y = mm 



E 



K 



i=l,...,Nb 
D=Di,...,Di, 



k=l 



k,k' = l 



(6) 



This generic definition encompasses all the spectral information [i index) from each detector {D index) 
and systematics parameterization through ak and Sf^. For Sf^ ~ 0, we recover the classical definition, 

xlo syst = E (^?)^' through 

Af = f^v;^ - TVf ) /[/f (7) 



where we assume that the simulated data event numbers, A'^ ^, arc uncorrelated between bins and detectors. 
In the absence of real data, they are computed for fixed given values of Am§i* and sin^(20i3)*. On the 
other hand , the theoretical model, relies on the searched sin^(20i3) value. We assume an uncorrelated 
weight error U[' which, in the absence of systematic uncertainty, is simply expressed as the statistical error: 

= 

Systematic uncertainties are included in the definition ([6|) through ak and Sf^ coefficients. The Sfj, 
coefficient represents the shift of the i"' bin of detector D spectrum due to a 1 cr variation in the fc'^ 
systematic uncertainty parameter ak- Following this definition, we introduce bin, detector and reactor 
correlations in the systematic errors through Sj^i. definitions whereas systematic parameter correlations are 
gathered in Wk^k' (we refer to the appendix for details). Eventually, some fully uncorrelated systematic 
uncertainties may be included through the UP definition, by adding quadratically all their effects together 
with the statistical uncertainty. We will use this property to include background shape uncertainties in 
our analysis as described in section [S] We refer to the appendix for the proper inclusion and definition of 
systematic coefficients 5*/^^. and ak inside the definition 

We define the sensitivity or sensitivity limit sin^ (26*13)11111 as the largest value of sin^ (26*13) which fits the 
true value sin^(20i3)* = at the chosen confidence level. We therefore determine the sin^(20i3) sensitivity 
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at 90 % C.L. of an experiment as the value of sin^(20i3) for which 

Ax'=x'(sm2(2^i3))-xLn = 2.71 



(8) 



4 Generic overview of systematic error inputs 

Systematic errors can be classified into three main categories: reactor, detector and data analysis induced 
uncertainties. In this section we provide a brief and generic description of the systematic uncertainties 
included in our modclization. Details concerning specific experimental cases are given in sections IHl [3 

The dominant reactor induced systematic error comes from our limited knowledge of the physical processes 
which produce electron antineutrinos in nuclear reactors. This leads to an overall normalization error on 
the production rate of 1.9 % jH]. Similarily we include a 2.0 % uncertainty on the antineutrino spectral 
shape [15] . with the conservative hypothesis that the energy bins are not correlated among themselves. 
Furthermore, even at a perfectly defined thermal power and with an absolute knowledge of the number of 
antineutrinos emitted for each fission, an underlying uncertainty remains since the nuclear energy released 
per fission is known to about 0.5 % [S]. In our model the last three uncertainties are taken to be fully 
correlated between the nuclear cores. 

We included another group of reactor induced systematics, taken to be uncorrelated between the reactor 
cores: the uncertainty on the determination of the thermal power of each nuclear core, within 0.6-3 % and 
the uncertainty on the isotopic composition of the nuclear core fuel elements, within 2-3 %. Also, finite 
size and solid angle effects, distances bias between reactors and detectors, as well as displacements of the 
neutrino production barycentcr might affect the fluxes at the near detector (s) if they are close to the power 
plants (below 200 m) up to a level of 0.1 %. 

We did not implemented the uncertainty coming from neutrinos produced in the spent fuel pools, often 
located within tens of meters from the nuclear cores since we could not gather all the relevant information 
for the different sites. We thus neglect this effect in our simulation, whithout any justification except in 
the case of Double Chooz, since a detailed evaluation has shown that this uncertainty does not affect its 
sensitivity [24]. We point out that in particular cases this additional neutrino source slightly affects the 
antineutrino spectrum, and is thus relevant for experiments aiming at high sensitivities, e.g. sin^(20i3) ~ 
0.01. We could also have introduced a specific error on the inverse-beta decay neutrino cross section of 
0.1 % mini]. However, being fully correlated between the detectors, the latter can be gathered into a global 
uncertainty, adding up to the overall neutrino rate knowledge. 

The basic principle of the multi-detector concept is the cancellation of the reactor induced systematics; 
additional contributions would not modify significantly the sensitivities of the forthcoming experiments. 
Let us now focus on the uncorrelated errors between detectors that could affect strongly the experimental 
results. 

The uncorrelated errors between detectors directly contribute to the relative normalization of the mea- 
sured antineutrino energy spectrum of each detector. One of the major improvement with respect to the 
CHOOZ experiment relies on the precise measurement of the number of free protons inside target volumes, 
proportional to the antineutrino rates. Experimentally the target mass will be determined at 0.2 %. An 
uncertainty on the fraction of free hydrogen per unit volume remains, at the level of 0.2-0.8 %, if different 
batches of liquid scintillator are used to fill the detectors of a given experiment. This complex case diserve a 
special treatment as described in section [7^ We did not include any error associated with the measurement 
of the live time of the experiments. 

The last set of systematics concerns the selection cuts applied to extract the antineutrino signal and reject 
the backgrounds. Neutrino events are identified as positrons followed in time by a single neutron captured on 
a gadolinium nucleus. New detector designs have been proposed in order to simplify the analysis, reducing 
the systematic errors while keeping high statistics and high detection efficiency. We accounted for three 
uncertainties for both positron and neutron associated to a candidate event: the possibility of missing the 
particle as it escapes the target (escape), the uncertainty related to the particle interactionqj (capture) and 

■^Note that it affects essentially the neutrons, since the Gd concentration might differ between the detectors if they are not 



4 



the identification cut based on the energy deposited in the detector active region (identification). Finally 
we take into account the uncertainty on the efficiency of the delay cut and an error on the neutron unicity 
of the event selections, whereas we do not consider any position vertex reconstruction. We provide in 
section l6TT] (Table [1]) the detailed systematic inputs necessary for our simulations. The uncertainty on these 
efficiencies have been treated, otherwise stated, as uncorrelatcd between detectors. Uncertainties induced 
by the background subtractions are discussed in sections [5] and 16.31 

5 Backgrounds: description and modelization 

In this section we briefly review the three main kinds of background for the next generation of reactor 
neutrino experiments: accidental coincidences, fast neutrons and the long-lived muon induced isotopes 
^Li/^He. We then describe our simplified background modelization. 

Naturally occurring radioactivity mostly creates accidental background, defined as a coincidence of a prompt 
energy deposition between 0.5 and 10 McV, followed by a delayed neutron-like event in the fiducial volume 
of the detector within a few hundredths of a millisecond. Selection of high purity materials for detector 
construction (scintillator, mineral oil, PMT's, steel, etc.) and passive shielding provide an efficient handle 
against this type of background. We assume that the accidental background rate can be measured in situ 
with a precision of 10 %. 

Cosmic ray muons will be the dominating trigger rate at the depth of all near detector sites. Even though the 
energy deposition corresponds to about 2 MeV per centimeter path length (providing a strong discrimination 
tool) they induce the main source of background. 

Muon induced production of the radioactive isotopes *He, ^Li and ^^Li can not be correlated to the primary 
muon interaction since their lifetimes are much longer than the characteristic time between two subsequent 
muon interactions. The characteristic signature of this last class of events consists in a four-fold coincidence 
(jU ^ n ^ /3 — *■ n). The initial muon interaction is followed by the capture of spallation neutrons within 
about 1 ms. The time scale of the decay of the considered isotopes is on the order of a few 100 ms, again 
followed by a neutron capture. This background mimics the 77e signal and is considered among the most 
serious difficulty to overcome for the next generation reactor neutrino experiments. In our simulation we 
will assume that it can be estimated to within 50 %. 

A further source of background are neutrons that arc produced in the surrounding rocks by radioactivity 
and in cosmic ray muon induced hadronic cascades. In the latter case, which is dominant at shallow depth, 
the primary cosmic ray muon may not penetrate the detector, being thus invisible. Fast neutrons may 
then enter the detector and create recoil protons and be captured by hydrogen or gadolinium nuclei after 
thermalization. Such a sequence can mimick a Fg event. In the case of Double Chooz far detector (depth of 
300 m.w.e.), muon induced neutron production can be fairly well estimated from the results of the CHOOZ 
experiment, since it was the dominating background monitored during reactor off periods ^ . We will assume 
that this background rate can be estimated within a factor of two. Figure [T] illustrates the background 
spectra that we implemented in our modelization. We used the CHOOZ reactor off data [9] to estimate the 
fast neutron background spectral shape, measured to be flat. We used a simple approximate exponential 
shape for accidental backgrounds. Finally we implemented the spectra of ^He and ^Li based on nuclear data 
information; we weighted them in the ratio 0.2/0.8, respectively. We checked that slight modifications of 
the background shapes do not change the results of our simulations for sin'^(26'i3) sensitivities above 0.01. 

Background rates are adjusted to match to the predictions at each detector location. Afterward they are 
subtracted from the total i?vis spectrum in each detector. These three backgrounds, i3„, have rates and 
shapes known at ago and (Tghp^^D , respectively. We take the conservative assumption that these shape 
uncertainties are fully uncorrelatcd between bins. Thus, the (Tshp.B'^ contributions may be directly included 
inside Uf (eq. 0). These generic uncorrelatcd errors will then take into account the statistical and the 
background shape subtraction uncertainties: 

filled with a single batch of liquid scintillator. 
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Figure 1: Energy spectrum of the backgrounds from spallation neutrons, accidentals, and cosmogenic 
20 % *He and 80 % ^Li. Each curve is normalized to unity. 
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(9) 



As the background rate uncertainties are correlated between bins, we treat them with additional pull terms, 
asn, in eq. with weights ago included in Sf'/, (see the appendix). 



6 Comparison of the current proposals 

Several sites are currently being considered for new reactor experiments to search for ^13: Angra dos 
Reis (Angra, Brazil), Chooz (Double Chooz. France, and possibly Triple Chooz), Daya Bay (Daya Bay, 
China), Kashiwazaki (KASKA, Japan) and Yonggwang (RENO, Korea). All these experiments may be 
classified in two generations. The first aims to probe the value of sin^(26'i3) till 0.02-03, and the second to 
track sin^(2^^i3) down to 0.01 (90 % C.L.). The first phase concerns Double Chooz, RENO, and possibly 
Daya Bay (with its phase I). This phase should end by 2013. Angra, Daya Bay (nominal setup with 
8 detectors), KASKA and possibly Triple Chooz are focusing on the second phase. For these second 
generation experiments, a significant R&D effort is required since the effective Gd-scintillator mass will be 
increased by, at least, one order of magnitude, and systematics, as well as backgrounds uncertainties, have 
to be further reduced with respect to the first phase experiments. In the following comparisons we will not 
include Angra, KASKA and Triple Chooz. We will only focus on Double Chooz, Daya Bay and RENO. 

In the next sections we first introduce a generic discussion on systematics, scintillator composition and 
backgrounds fsections 16.11 to [H75|l . We will then shortly describe each setup and compute the associated 
baseline sensitivity fsections 17.11 to [775]) . We then perform a comparative analysis of the setups in section [S] 
to highlight advantages and drawbacks of each setup. 



6.1 Detailed systematics review 

The two Double Chooz and Daya Bay proposals take a careful inventory of systematics, compiled 
in Table [T] for comparison. Double Chooz and Daya Bay estimates for the case of no additional R&D are 
at hand. We decided to strictly use the systematic errors (Table [T]) and background values quoted by the 
collaborations [24l ESI [26] . However, we could not find any detailed background estimate for the RENO and 
Daya Bay Mid site setups. In the latter cases wc use a simple model to estimate the background subtraction 
uncertainties from the scaling of the Double Chooz far detector (see section 1575)) . 
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Error Description 


CHOOZ 


Double Chooz 




Day a Bay 














No R&D 


R&D 




Absolute 


Absolute 


Relative 


Absolute 


Relative 


Relative 


Reactor 








Production cross section 


1.90 % 


1.90 % 




1.90 % 






Core powers 


0.70 % 


2.00 % 




2.00 % 






Energy per fission 


0.60 % 


0.50 % 




0.50 % 






Solid angle/Bary. displct. 






0.07 % 




0.08 % 


0.08 % 


Detector 








Detection cross section 


0.30 % 


0.10 % 




0.10 % 






Target mass 


0.30 % 


0.20 % 


0.20 % 


0.20 % 


0.20 % 


0.02 % 


Fiducial volume 


0.20 % 












Target free H fraction 


0.80 % 


0.50 % 




? 


0.20 % 


0.10 % 


Dead time (electronics) 


0.25 % 












Analysis (paticle id.) 








6+ escape (D) 


0.10 % 












e+ capture (C) 














e+ identification cut (E) 


0.80 % 


0.10 % 


0.10 % 








n escape (D) 


0.10 % 












n capture (% Gd) (C) 


0.85 % 


0.30 % 


0.30 % 


0.10 % 


0.10 % 


0.10 % 


n identification cut (E) 


0.40 % 


0.20 % 


0.20 % 


0.20 % 


0.20 % 


0.10 % 


Fe time cut (T) 


0.40 % 


0.10 % 


0.10 % 


0.10 % 


0.10 % 


0.03 % 


Fe distance cut (D) 


0.30 % 












unicity (n multiplicity) 


0.50 % 








0.05 % 


0.05 % 


Total 


2.72 % 


2.88 % 


0.44 % 


2.82 % 


0.39 % 


0.20 % 



Tabic 1: Breakdown of the systematic errors included in the computation of the sensitivity of Dou- 
ble Chooz [21] and Daya Bay [25]. Since no breakdown of the RENO systematic errors has been pub- 
lished we use the same systematic error budget as for Double Chooz. Double Chooz and Daya Bay relative 
systematics are almost comparable, the only main difference coming from the determination of the gadolin- 
ium concentration and the free proton fraction inside the target volume. The absolute determination of 
the free proton fraction can have some impact in Daya Bay since multiple batches will be used to fill all 
the 8 detectors (see section [T^ for details). Nevertheless, there is no pubhshed value in the Daya Bay 
proposal 
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Through all the available publications [211 US [HI [55] , the differences between the systematics are found 
only for the relative normalization of the dectetor ((Trei) and the subtraction uncertainties on background 
rates (ctbd). From section [4| and Table [H we thus group systematics in two categories to be used in our 
analysis eq. ([6|): 

1. generic systematics common to all the experiments (Tabic [2]): 

~ f abs , the theoretical uncertainty on reactor antincutrino spectrum prediction. We call it also the 
absolute normalization of event rates (since common to all the detectors), extracted from Table [T] 
without power uncertainty contributions, dabs is at the level of 2 %; 

- CTshp, the theoretical reactor spectrum shape uncertainty, at the level of 2 %; 

- (7sci, iTg^J, the absolute and relative energy scale uncertainties, roughly assessed at the level of 
0.5 % each; 

~ CTpwr, the reactor thermal power uncertainty, at the level of 2.0 %; 

- (Temp, the reactor core specific fuel composition uncertainty which is roughly at the level of the 
power uncertainties (2-3 %) on each fuel clement. 



Cabs 


Cshp 


_abs 


rol 


Cpwr 




2.0 % 


2.0 % 


0.5 % 


0.5 % 


2.0 % 


2-3 % 



Table 2: Generic systematic uncertainties as included in the analysis. For more details about the 
correlations between the systematics, we refer to the appendix, and more particularly to Table [T5l 



2. specific systematics: 

- fTrci, the relative normalization of event rates between all the detectors. This uncertainty is 
uncorrclated between detectors; 

- CT^-D, the background subtraction unceratinties, described in section [6.31 

6.2 Impact of the scintillator composition 

In this section we stress the impact of the scintillator composition on the sensitivity of reactor neutrino 
experiments. All current projects will use a gadolinium doped liquid scintillator to enhance the neutron 
capture. The long-term stability of this scintillator is among the most difficult experimental challenges, 
since a degradation of the scintillator transparency would induce large systematic uncertainties. In the 
following we consider a stable scintillator for all experiments. The choice of the scintillating base has 
some importance since it defines the free proton number per unit volume, the rate and proton recoil 
background rate. Similarly the ^^C number per volume drives the long-lived muon induced isotopes in the 
target scintillator. 

Different bases can be used as neutrino target scintillator, mixture of dodecane (DOD), pseudocumene 
(PC) or phenybcylylethane (PXE), or linear alkylbenzene (LAB), as described in details in Table [H If we 
consider the Double Chooz scintillator (80 % DOD -|- 20 % PXE) as our reference, a pure LAB scintillator 
contains 4.9 % less free proton per volume, and 9.5 % more carbon atoms. In the following we will assume 
the Daya Bay [2S] and RENO [5S] experiments will use pure LAB as the target scintillator, and we will 
renormalize the neutrino and the backgound rates accordingly in section 17.21 and 17.31 

6.3 Estimation of the ^u-induced backgrounds 

Our estimates of the backgrounds induced by cosmic muons is based on the Double Chooz proposal [21] . A 
modification of the MUSIC code [5D| was used to compute the muon rates and energy spectra by propagating 
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Liquid 


Formula 


density 


102« H/m3 


1028 C/m3 


Dodecane (DOD) 


\ 1 1 1 1 /I? 


0.753 


6.93 


3.20 


nQf^n fi riPiiTnono i Pi! 1 


C9H12 


u .00 


5 30 


3 Q7 








5 08 

(J .WO 


4.52 


90 % DOD+10 % PC 


m 1 1 1 vo i n vol 


0.77 


6.77 


3.32 


80 % DOD+20 % PXE 


mixture in vol. 


0.80 


6.56 


3.46 


linear alkylbenzene (LAB) 


C16H30 


0.86 


6.24 


3.79 



Table 3: Impact of the scintillator composition on the target free proton number driving the rate 
as well as the recoil proton background, and on the carbon composition which drives the long-lived muon 
induced isotopes produced in the detector. In the following wc consider that Double Chooz is using modules 
containing 8.26 tons of 80 % dodecane +20 % PXE (in volume) based target scintillator, and that Daya Bay 
and RENO are both using 20 ton modules of a LAB based target scintillator. 



surface muons through rock. The site topographies have been included, according to a digitized map of the 
Chooz hill profile [12] in the case of the far site (300 m.w.e. overburden), and according to a flat topography 
in the case of the near site (80 m.w.e. overburden). For the far site this full Monte-Carlo simulation predicts 
a muon flux <I>^=0. 612 ± 0.007 m~2s~^, slightly higher than the approximate measured value quoted in [9]. 
The mean muon energy computed according to this method is (E^) — 61 GeV. For the near site, we get 
(f>^ = 5.9 m~^s~^, and (i?^) = 22 GeV. A similar detailed computation was performed for the three sites 
of the Daya Bay collaboration [15] . 

Thus, for both the cases of Double Chooz and Daya Bay we use only the values of the muon flux and mean 
energy computed by the collaborations. However, we could not find any published data for the case of 
RENO. We then use the underground muon fluxes and mean energies calculated analytically following [IT] . 
In order to justify this approximation. Table H] reports the muon flux and mean energy computed by the 
Double Chooz and the Daya Bay collaborations, from 300 m.w.e. to 923 m.w.e., as well as the analytical 
computations according to [21] . We found a reasonable agreement which bears out the use of the analytical 
model [H] for the case of the RENO sites (depths of 255 m.w.e. and 675 m.w.e.). Nevertheless, note 
that the analytical simulation assumes a flat topography. It is worth noting, however, that the mean muon 
energy predicted by the analytical computation is systematically ~ 25 % lower in the depth range of interest. 
Therefore we arbitrarily renormalize the analytical calculation by 25 % to estimate the mean muon energy 
at the RENO sites. We also note a large discrepancy between the detailed computation and the analytical 
model of [H] for the Double Chooz near site, probably due to its very shallow depth. 

Cross sections of muon induced isotope production on liquid scintillator targets (^^C) have been measured 
by the NA54 experiment at the CERN SPS muon beam at 100 GeV and 190 GeV muon energies [IHj. The 
energy dependence was found to scale as atot{E^) (x (i?^)" with a — 0.73 ± 0.10 averaged over the various 
isotopes produced. We consider in the following that both long-lived muon induced isotopes and muon 
induced fast neutron backgrounds scale as 'I'^ x (Ef^)", our reference being taken at the Chooz far site (full 
Monte-Carlo simulation). We define the depth scaling factor as 

DSF = ($, {E,r)/i^, {EX) Double Chooz far ? 

(10) 

which is illustrated in the last column of Table|3|for various detector sites. Daily background rates computed 
for Daya Bay, Double Chooz, and RENO detectors at the different sites are then given in Table jSj 

Two cases are considered and compared: the background rates taken from the literature when available, 
and the background rates extrapolated from the background computed for the Double Chooz far detector, 
scaled with the target mass, the scintillator free proton and carbon numbers, as well as the depth scaling 
factor (DSF). We note here the good agreement between the Daya Bay cosmogcnic induced backgrounds 
(^Li/^He and fast neutrons) estimated from our Double Chooz extended model and the original estimates 
of the Daya Bay collaboration. This bears out the use of our model for the RENO and Daya Bay mid 
site configurations. For the case of the ^Li/^Hc background, we understand well this agreement since 
the background rate mainly depends on the mass of the neutrino target region. The agreement is more 
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Site 


deptn (m.w.e.j, 


Detailed simulation 


Analytical 


model 


DSF 




topography 
















m"2 s^i 


GeV 


m ^ s ^ 


GeV 




DC near 


80, flat 


5.9 


22 


9.9 


17 


6.80 


RENO near 


230, hill 






1.2 


40 


1.57 


DB near 1 


255, hill 


1.2 


55 


0.9 


44 


1.32 


DB near 2 


291, hill 


0.73 


60 


0.72 


49 


1.06 


DC far 


300, hill 


0.61 


61 


0.67 


50 


1 


DB mid 


541, hill 


0.17 


97 


0.15 


71 


0.32 


RENO far 


675, hill 






0.084 


94 


0.20 


DB far 


923, hill 


0.04 


138 


0.035 


118 


0.10 



Table 4: Muon flux <i>^ and mean energy {E^j) for the underground site of the reactor neutrino experiments. 
We compare the values obtained from a full Monte-Carlo simulation for Double Chooz and Daya Bay to the 
analytical model of [21]. We use the latter model for RENO. The depth scaling factor (DSF) is defined by 
the product x {E^)°')/ x (iS^^)") Double Chooz far, the Doublc Chooz far site is taken as the reference. 
Backgrounds induced by cosmic muons are scaled according to this factor. 



Detector 


Accidental (d ^) 


/^t-induced fast-n (d ^) 


/^-induced 


9Li/«He (d-i) 


Site 


Original 


DC ext. 


Original 


DC ext. 


Original 


DC ext. 


Double Chooz near 




13.60 ± 1.36 




1.36 ±1.36 




9.52 ±4.76 


RENO near 




7.10 ±0.71 




0.68 ± 0.68 




5.40 ± 2.70 


Daya Bay DB 


1.86 ±0.19 


5.98 ±0.60 


0.50 ±0.50 


0.57 ±0.57 


3.7 ± 1.85 


4.55 ±2.27 


Daya Bay LA 


1.52 ±0.15 


4.76 ± 0.48 


0.35 ±0.35 


0.45 ± 0.45 


2.5 ± 1.25 


3.63 ±1.81 


Double Chooz far 


2.00 ±0.20 




0.20 ±0.20 




1.40 ± 0.70 




Daya Bay mid 




1.45 ±0.14 




0.14 ±0.14 




1.10 ±0.57 


RENO far 




0.90 ±0.09 




0.09 ± 0.09 




0.69 ±0.35 


Daya Bay far 


0.12 ±0.01 


0.44 ± 0.04 


0.03 ± 0.03 


0.04 ± 0.04 


0.26 ±0.13 


0.33 ±0.17 



Table 5: Daily background rates computed for Daya Bay, Double Chooz, and RENO detectors at the 
different sites. We consider the three main background sources: accidental events, /i-induced fast neutrons, 
and ^-induced ^Li/^He. The columns labelled "Original" quote the background rates taken from the 
literature when available. The columns labelled "DC ext." (for extended) quote the background value 
extrapolated from the background computed for the Double Chooz far detector, scaled with the detector 
target mass, the scintillator free proton and carbon numbers, as well as the depth scaling factor (DSF). 
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Background subtraction systematic errors (%) 



Detector 
Site 


Accidental 
Original DC ext. 


/i-induced fast-n 
Original DC ext. 


/^t-induced '■*Li/**He 
Original DC ext. 


Double Chooz near 




0.123 




0.123 




0.043 


rLH/INO near 


— 


0.019 


— 


0.019 


— 


0.074 


Uaya oay Dt> 


0.020 


0.064 


0.054 


0.061 


0.199 


0.245 


Uaya rJay LA 


0.020 


0.063 


0.046 


0.060 


0.164 


0.239 


Double Chooz far 


0.292 




0.292 




1.020 




Daya Bay mid 




0.120 




0.115 




0.458 


RENO far 




0.100 




0.095 




0.382 


Daya Bay far 


0.010 


0.036 


0.025 


0.035 


0.108 


0.138 



Table 6: Background subtraction systematic errors (in percent) computed for Daya Bay, Double Chooz, and 
RENO detectors at the different sites. We consider the three main background sources: accidental events, 
/x-induced fast neutrons, and ^-induced ^Li/^He. The columns labelled "Original" quote the systematic 
errors taken from the literature when available. The columns labelled "DC ext." (for extended) quote 
the systematic errors extrapolated from the Double Chooz far detector, taking into account the estimated 
detector signal to background ratio as well as the background rate uncertainty. 



surprising, however, for the case of the fast neutron background, since the size of the liquid shielding around 
the detector active area is rather different between the Double Chooz and Daya Bay detectors (the RENO 
detector design is very close to the Double Chooz case). In addition, further detector differences such as 
the thickness of the buffer oil shielding the inner target region, as well as the different mechanical structure 
explain the discrepancy between the Daya Bay computation and the DC extended model for the case of 
the accidental background. Nevertheless, we found out that these differences influence only weakly the 
sensitivity computed for the three experiments since the accidental background energy spectrum is different 
enough from the expected oscillation signal, and it is supposed to be known with a precision of 10 %. 

In a similar way. Table [S] gives the background subtraction systematic errors (in percent) computed for 
Daya Bay, Double Chooz, and RENO detectors at the different sites, taking into account the estimated 
detector signal to noise ratios as well as the background rate uncertainties. 

7 Reactor experiments baseline sensitivity 
7.1 Double Chooz 

The Double Chooz collaboration is composed of institutes from Brazil, France, Germany, Japan, Russia, 
Spain, United Kingdom, and the United States. The experimental site is located close to the twin reactor 
cores of the Chooz nuclear power station (two PWb[1 producing 8.5 GWth), operated by the French company 
Electricite de France (EDF). The two, almost identical, detectors will contain a 8.3 ton fiducial volume of 
liquid scintillator (density of 0.8) doped with 0.1 g/1 of gadolinium (Gd). The far detector will be installed 
in the existing laboratory, 1.05 km from the cores barycenter, shielded by 300 m.w.e. of rock. This detector 
should be operating alone for 1.5-2 years (DC Phase I), starting data taking by the end of 2008. The 
second detector will be installed in the meantime about 280 m from the nuclear core barycenter, at the 
bottom of a 40 m shaft (80 m.w.e.) to be excavated. Distances between detectors and nuclear cores as well 
as site overburdens are given in Table [71 Since there are no more than two NPP cores, it is still possible to 
install the near detector at a suitable position where the ratio of reactor Ve fluxes from each core is the same 
as for the far detector (the iso-ratio curve is plotted on figure [5]). This allows reactor relative uncertainty 
cancellations (NPP core compositions). This detector should be operational by 2010, and will take data 

^Pressurized Water Reactor 
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Detector 


near 


far 


Distance from West reactor (m) 


290.7 


1114.6±0.1 


Distance from East reactor (m) 


260.3 


997.9 ±0.1 


Detector Efficiency 


80 % 


80 % 


Dead Time 


25 % 


2.5 % 


Rate without efficiency (d~^) 


977 


66 


Rate with detector efficiency (d~^) 


782 


53 


Integrated rate (y""^) 


1.67 10^ 


1.4810-* 



Table 7: Double Chooz antineutrinos rate expected in the near and far detectors, with and without reactor 
and detector efficiencies. The integrated rate in the last line includes detector efficiency, dead time, and 
reactor off periods averaged over a year. The averaged reactor global load factor is estimated at 79 % [3D] . 



\ 



\ 



280 m 



V 4 




•f^ Reactor 
O Detector 



1 ,050 m 



Iso-ratio curve 



Figure 2: Double Chooz experiment site configuration. We show also on this figure the far flux iso-ratio 
line. Another detector located on this particular curve will receive the same reactor flux ratio as for the far 
detector: 44.5 % from West and 55.5 % from East reactor. The near detector of Double Chooz is foreseen 
to be placed on this line. 
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for three years (DC Phase II). Other details concerning the experiment may be found in the collaboration 
proposals 

For the Double Chooz phase I (DC Phase I) analysis, we used systematics of Table [U setting (Troi and 
"'sci to since only one detector will be present. For a data taking period of 1.5 years, the sensitivity 
is sin'^ (26*13)11111 = 0.0544. The second phase (DC Phase II) will then start and both detectors will take 
data for 3 years as scheduled in the proposal. The full experiment will then achieve a final sensitivity 
sin^(20i3)ii,ii = 0.0278, assuming systematics from Table [2] and (Trci = 0.6 % j24]. The sensitivity worsens 
for Arri^i < 2.5 10~^ eV^, due to the close distance of the far detectoi0. If we take the lower bound [121 113] 
on Arrig]^ (2.0 10~^ eV^), Double Chooz will lose 30 % in sensitivity whereas for the upper bound [T^ [T^ 
on Amo, (3.0 10-3 eV^), Double Chooz will gain 15 % in sensitivity. 

7.2 Daya Bay 

Daya Bay is an experiment proposed by institutes from China, the United States, and Russia. Daya Bay |25| 
will be located in the Guang-Dong Province, on the site of the Daya Bay nuclear power station. The site is 
made up of two pairs of twin reactors, Daya Bay (DB) and Ling Ao I (LA I). An additional pair of reactors. 
Ling Ao II (LA II), is currently under construction and should be operational by 2010-2011 [H]. Each 
core has a thermal power of 2.9 GW [30]. In the full installation setup 3.3 km of tunnel and 3 detector halls 
have to be excavated, in order to accommodate 8 detector modules [25]. Each module contains an effective 
volume of 20 tons of Gd- loaded LAB liquid scintillator (Table [3]). Distances between detectors and nuclear 
cores as well as site overburdens are given in Table [S] This site yields a rather complex signal composition 



Detector 


near DB 


near LA 


mid 


far 


Distance from DB 1 (m) 


350 


1,356 


1,153 


1,970 


Distance from DB 2 (m) 


381 


1,331 


1,161 


2,000 


Distance from LA I 1 (m) 


942 


492 


783 


1,619 


Distance from LA I 2 (m) 


1,030 


475 


818 


1,623 


Distance from LA II 1 (m) 


1,378 


500 


968 


1,602 


Distance from LA II 2 (m) 


1,463 


555 


1,029 


1,624 


Detector eff. 


80 % 


80 % 


80 % 


80 % 


Dead Time 


7.2 % 


4.3 % 


1 % 


0.2 % 


Rate without eff. (d"^) 


1,938 


1,813 


494 


430 


Rate with detector eff. (d^^) 


1,550 


1,450 


395 


344 


Integrated rate (y"^) 


4.10 10^ 


3.95 10^ 


1.11 10^ 


9.7710'* 



Tabic 8: Daya Bay antineutrino rates expected in the near and far detectors, with and without reactor and 
detector efficiencies. The integrated rate in the last line includes detector efficiency, dead time, and reactor 
off periods averaged over a year. We assumed that LA II NPP will be operating for the time the far site 
will be fully installed, but for the Mid site installation we assumed that LA II will be off. 

in each detector coming from up to 6 different NPP cores, as shown in Table [51 According to the Day Bay 
proposal, our estimate of the sensitivity is sin^(20i3) = 0.0085 very close to the Daya Bay quoted value 
(sin^(20i3) = 0.008). However, in the Daya Bay proposal, the uncertainties on the reactor fuel composition, 
fcmp, as well as the energy scale associated uncertainties a^^f and (Tg°} are neglected. Taking into account 
these systematics (Table [2]) and the quoted value of ctioi = 0.39 % with no R&D [25], our estimate of the 
Daya Bay final sensitivity is then sin^(20i3) = 0.009, with the full installation (DB Phase 11, Figure [3]) 
after 3 years of data taking. We draw the attention of the reader on the point that these computations are 
based on the assumption that cFi-d ~ 0.39 % is fully uncorrelated between all the detectors, and in particular 

''a bit too close to the NPP to get the maximum amount of information from first minimum of oscillation over the reactor 
spectrum distortion 
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DB LAI LA2 



near 1 (DB) 
near 2 (LA) 
Mid (LA2 OFF) 
Mid (LA2 ON) 
far 



83.1 % 
6.5 % 
32.3 % 
22.5 % 
24.9 % 



11.4 % 

50.6 % 

67.7 % 
47.1 % 
37.4 % 



5.5 % 
42.8 % 
0.0 % 
30.4 % 
37.7 % 



Table 9: Daya Bay rate contributions from each NPP set (2 cores by set) while assuming, except if otherwise 
noticed (3''^ line), the new Ling Ao II NPP is operating at full power. 



between detectors on a same site. Nevertheless this hypothesis is not guarenteed. We discuss this point as 
well as the current filling scenario (multi-batches) just below. 



Figure 3: Daya Bay installation phases and site configuration. On the left: phase I, with 2 x (2 x 20 t) 
detectors, Daya Bay and Ling Ao I power plants are operating. On the right: phase II with 2 x (2 x 20 t) 
near sites, and 4 x 20 t at far site, all three power plants are operational. 

A preliminary fast phase (DB Phase I) is proposed by the Daya Bay collaboration. This phase includes only 
4 detectors, 2 of them located at the DB near site and the other 2 at the mid site (see Figure Taking 
systematics from Table [5] and a^d = 0.39 % fully uncorrelated between detectors, we get a sensitivity after 
1 year of data taking of sin^ (26*13) = 0.040. if LA II NPP is off. If DB Phase I starts after 2010, with LA II 
operational, we get a sensitivity of sin^ (26113) = 0.038 still after 1 year of data taking. The reason for a 
better sensitivity in the second scenario is that the mid detectors get 44 % more events (Table [51 [5] and 
Figure [3]), with a larger oscillation baseline with respect to LA I. 

Daya Bay site correlation 

The main concept of the Daya Bay experiment is based on the multi-inter-calibration of detectors. Since 
many detectors are installed on a same site, the total uncorrelated uncertainty of a site is decreased by 
a factor of l/VWf compared to the single detector uncorrelated uncertainty, where is the number of 
detectors on a given site. In the Daya Bay proposal, the full relative uncertainty {a^ei = 0.39 %) is assumed 
to be uncorrelated between all the detectors. However, the fraction of correlated and non-correlated error 




DB 



DB 
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between detectors of a same site is not trivial. It relies on many experimental assumptions on uncertainty 
correlations. The absence of correlations between detectors on a same site implies there will be no detector- 
to-detector correction applied. If detector responses happen to be different, any data correction from 
detector to detector on a same site would yield correlations between detector uncertainties. 

For DB Phase II, if we assume that CTroi is fully uncorrelated between the detectors, we get a sensitivity of 
sin^(20i3) = 0.009. On the contrary, if we assume that (t,.o1 is fully uncorrelated between detectors in different 
sites, but fully correlated between detectors on a same site, we get a sensitivity of sin'^(2^^i3) = 0.012. The 
real sensitivity should lay between these two extremes. 

For DB Phase I, if we assume CTi-cI is fully correlated on a same site and fully uncorrelated on distant sites, 
we get sin2(26li3) = 0.041 if LA II is off, and sin2(26li3) = 0.038 if LA II is on. We conlude from these results 
that DB Mid sensitivity does not strongly depend on the correlations in CTi-cI. More generally, DB Mid setup 
weakly depends on a^ci- This latter point will be discussed in section [H together with a full description of 
the impact of each systematic on the forseen sensitivity. 



Daya Bay and the filling procedure 

A large amount of Gd-loadcd liquid scintillator will have to be produced, stored and filled into the detectors 
(8x20 tons, which is 10 times more than in Double Chooz). Due to the large number of detectors and 
the large amount of liquid to manage, the Daya Bay collaboration plans to fill detectors with four different 
Gd-doped liquid scintillator batches. A single batch will be used to fill detectors by pairs (Figured]). 
The best installation scenario (which is the one chosen by the Daya Bay collaboration) is then to move 
one of the filled detector to a near site, the other one to the far site (scheme 1 of Figure |3] and Table [TO]). 
Another batch of Gd-loaded liquid scintillator will be used for the next pair of detectors and so on. With 
the adopted filling procedure [25j , extra systematic uncertainties on the hydrogen content between different 
batches have to be included. 



# Pair 1 
■ Pair 2 
A Pair 3 
■uT Pair 4 



Far 


DB 


LA 



Scheme 1 



Scheme 4 





•■ 


AA 



Scheme 2 





•A 





Scheme 5 





■A 


■A 



Scheme 3 





■■ 


At? 



Scheme 6 





mm 


AA 



Figure 4: Possible installation of detector pairs in the Daya Bay experiment according to the adopted filling 
procedure [IS]. Due to the large number of detectors, and the large amount of liquid to be managed, 
the Daya Bay collaboration plans to fill detectors by pairs with four different Gd-doped liquid scintillator 
batches. The forseen detector installation scenario j^S] corresponds to scheme 1. We illustrate here other 
installation possibilities (schemes 2 to 6). 



Although the relative uncertainty on the free proton content within a same batch could be kept at a very 
low level (0.2 % in the Daya Bay proposal [5Sj, negligcable in the Double Chooz proposal [M]), it is not 
necessarily true in the case of different batches, in which the chemical composition may slightly change. 
The free proton content between different batches relies then on the measurement of this quantity. In the 
CHOOZ experiment [S], the free proton fraction inside the Gd-loadcd liquid scintillator was known to 0.8 %. 
In Double Chooz, this uncertainty is assessed at the level of 0.5 %. Since there is no published value on the 
absolute determination of this quantity in the Daya Bay proposal, we assume here, as in Table [1] a 0.5 % 
uncertainty between different batches. The filling systematic coefficients are explained in Table [10] (refer to 
the appendix and to Table [15] for full description) . The uncertainty on the free proton fraction of a single 
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Scheme 1 
Scheme 2 
Scheme 3 
Scheme 4 
Scheme 5 
Scheme 6 



0.008 



0.01 



Single batch 



J \ L 

0.012 



Four batches 



\ 

0.014 



0.016 



sin (29^3) 



Figure 5: Sensitivity on sin^ (26*13) at 90 % C.L. after 3 years of data taking for the 6 different installation 
schemes illustrated on FigureSl with cTpair = 0.5 %. Left bounds are computed assuming uncorrelated errors 
between all the detectors and right bounds are for the assumption of full correlation of CTroi between detectors 
of a same site. The real sensitivity should be somewhere in between these two bounds. For comparison, 
we also show on this graph the computations for the single batch hypothesis where we take a^oi = 0.39 %. 
Note that obviously the first installation scheme provides the best sensitivity. In these results we do not 
include any detector swapping scenario. 



Error type k S^^ X 

Filling {N, = 5Na + iVb + 5N, 
of batch 1 A^s +1 
of batch 2 Ns +2 
of batch 3 A^s +3 
of batch 4 Ns +4 

Table 10: Daya Bay, phase II (8 detectors), specific filling systematic parameters table in addition to 
standard one (see the appendix for details). Here we adopt a fully uncorrelated uncertainty of the different 
batches, Cpair = 0.5 %, but fully correlated in a same batch. 



2) 


'n^' - 




apair,l 


•^pair 






'^pair,2 


•^pair 






'^paiijS 


•^pair 






C^paii-,4 
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batch is taken to cTpair = 0.5 %. This uncertainty is taken to be fully correlated when detectors are filled 
with the same batch and fully uncorrelated otherwise. 

According to this filling procedure, the final sensitivity of DB Phase II would be sin^(20i3) = 0.0093 
instead of 0.0089 with the initial installation scenario (scheme 1 of Figure 2]). In all the other configuration 
schemes (2-6), which allows comparing on a same site at least two detectors filled with the same batch, 
the sensitivity is more largely weakened (Figured]), to an extent depending on (Tpair- 

Note that we did not include any detector swapping option in previous conclusions. In the Day Bay proposal, 
the retained installation scenario is the first scheme of Figure 2) The baseline swapping option is then the 
permutation of two detectors filled with the same batch, in 4 steps, 1 step per batch. On the one hand, 
the drawback of such a swapping scenario is that two detectors filled with the same batch will never be 
directly compared. On the other hand, configuration schemes 2-6 allow detector intercalibration within a 
pair. However, it should be noticed that the time spent in any configuration different from scheme 1 may 
decrease the combined final sin^ (26*13) sensitivity (Figure O. 

7.3 RENO 

The RENO experiment [25] will be located close to the Yonggwang nuclear power plant in Korea, about 
400 km south of Seoul. The power plant is a complex of six PWR reactors, each of them producing a thermal 
power of 2.74 GW |3D]. The Yonggwong power station is ranked number 4 in the world, with a total thermal 
power of 16.4 GW. Its power rating is often cited as an advantage of the RENO experiment. These six 
reactors are equally distributed on a straight segment spanning over 1.5 km. The average cumulative 
operating factors for the reactors are all above 80 %. Figure El shows the foreseen layout of the experimental 
site. The near and far detectors will be located 150 m and 1,500 m away from the center of the reactor 




Figure 6: RENO experiment site configuration. 

row (Tabic [TT]) . and will be shielded by a 88 m hill (230 m.w.e.) and a 260 m "mountain" (675 m.w.e.) 
respectively. Two neutrino laboratories have to be excavated and equipped in order to host the detectors. 
They will be located at the bottom of two tunnels having a length of 100 m and 600 m for near and far 
detector, respectively. In this configuration, the flux contribution from Ri to Rq ranges from 3 % to 39 % 
for the near detector, and from 15 % to 18 % for the far detector (Table [12]). 

Assuming the sytematics of Table[T]and, lacking of published data, fixing (Trci = 0.6 % as for Double Chooz, 
we obtain a final sensitivity after 3 years of data taking of sin^(26'i3) — 0.021. This is calculated for 
2 X 20 t detectors. However, a recent talk [27], quotes a different detector size, with 15 t of Gd-dopcd liquid 
scintillator. In that case, the final sensitivity after 3 years, would be sin^(20i3) ~ 0.023. 

As seen from Table [T^] the near detector monitors mainly the two central NPP cores, and the experiment 
sensitivity is quite affected by the reactor power uncertainties. If we compute the sensitivity with only the 
2 central NPP cores, we get even better results compared to the full NPP configuration. This means that 
four of the six cores are useless for the sin^(26'i3) measurement. This can be understood by the fact that the 
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Detector 


near 


far 


Til ci" ct n nc^ TTT\Yn T? 1 I m 1 




1 fi77 
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474 




1 )lQl"fl'nr'0 TTrMTl rn I TTl i 
J_^iO LdlH^tj 11 Ulll J_ VjO I ill J 


212 


1 'i07 


J_^lD LdlH^Lj 11 Ulll _L VjT: I 111 J 


212 


1 'i07 


llicj'l'finr'p irmn R (ml 

J— /lO IjCllUjC 11 Ulll J.\jfJ I 111 I 


474 




Tit cf" ct n r'f^ TTr\m T? h ( m 1 
l-'loLcllUjO llUill l.\U 
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ou /o 


ou /o 


Tlcici/n T^im^i 
J-zCdU. i llllC 




u.o /o 


T?^itp witlimit pffipipTiPv TH""'^^ 


2859 


121 


Rate with detector efficiency (d~^) 


2287 


97 


Integrated rate (y~"^) 


6.20 10^ 


2.8210'' 



Table 11: RENO antincutrino rates expected in the near and far detectors, with and without reactor and 
detector efficiencies. The integrated rate in the last line includes detector efficiency, dead time, and reactor 
off periods averaged over a year. 



Ri 


i?2 


i?3 


i?4 


i?5 


i?6 


near 3.0 % 
far 14.8 % 


7.8 % 
16.9 % 


39.2 % 

18.3 % 


39.2 % 

18.3 % 


7.8 % 
16.9 % 


3.0 % 
14.8 % 



Table 12: RENO rate contributions from each reactor core. 



sensitivity gained by statistics is compensated by a loss of information on the rates and energy spectra. 
Thus the appeal of this site is diminished. The sensitivity is equivalent to a 5.8 GWth NPP reactor neutrino 
experiment for cTpwr = 2.0 %. 

The RENO collaboration considers using 3 small very near detectors (200 - 300 kg) to monitor sub-groups 
of cores of the NPP. However, taking into account current knowledge on reactor spectra (dabs = 2.0 %, 
(Tshp = 2.0 % [H]), even dedicated detectors with 10^ events will not improve the thermal power knowledge 
below ~2 %. Thus, the overall sin^(20i3) sensitivity will not improve. 

8 Discussion 

We develop here a two step comparison of the experiments described before: Double Chooz phase I (single 
detector, DC I), phase II (both detectors, DC 11), Daya Bay phase I {DB Mid) and phase II {DB Full) and, 
RENO {RN). The first elements of comparison are based on a single core equivalent approach. Although 
purely hypothetical, this analysis provides a lot of information on the impact of the layout of the site 
(locations of NPP reactor cores and detectors) . The second approach, giving far more information on the 
impact of systematics on each experiment, is based on the pulls- approach, presented in section[3]and detailed 
in the Appendix. Note that in the following discussion we adopt the approximation that all the NPP cores 
operate with the same average efficiency. However, this assumption is not guaranteed, and, especially for 
NPP with many reactors such as RENO and Daya Bay experiments, running time and procedure of each 
NPP core have to be taken into account in the final analysis. For the Double Chooz experiment, which 
places the near detector on the flux iso-ratio line of the far detector, the full running operation time and 
procedure of each core is not needed to perform the final analysis. 
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8.1 The single core equivalent approach 



First of all we may simplify each experiment to its roughest single core equivalent (SCE) with total matching 
power P = J2r=i ^^ with P"^ the power of the r*'^ reactor and the number of available reactors on site. 
In this case, we compute the sin^ (26*13) sensitivity for each experiment (see Table [T3|) for their baseline 
option [551 nil US] but also for a single core equivalent, and the averaged near and far detector locations 
computed with 



1 



^1/2 



This average distance, L, yields the same event rates associated to near and far detectors of each experiment 
except for the Daya Bay Phase II experiment. For this particular case, we have to compute this average 
distance in two steps for the equivalent near detector. In this special case, since there are two near sites, 
we compute the L for the DB near site, L^b, and for the LA near site, Lla, and then compute the overall 
single near detector equivalent distance as 

Ln ^ [ ::Ldb'^ + i^Lla'] ■ (12) 



Following this SCE simplification, numerical values of sin^(20i3) sensitivity are gathered in Table [01 As 





DC 


DB Mid 
LA II OFF LA II ON 


DB Full 


RN 


jC (in r.n.u.) 


17 


67 


101 


303 


71 


Lfar (in m) 


1,051 


931 


951 


1,716 


1579 


Lncar (in m) 


274 


484 


576 


441^ 


325 


sin^(26'i3)iim 


0.0278 
0.0274 


0.0410 
0.0274 


0.0381 
0.0289 


0.0110 
0.0105 


0.0213 
0.0176 



Table 13: In this table we provide the luminosity (expressed in r.n.u., see eq. ©), sin2(20i3)ii„, at 90 % C.L. 
of Double Chooz {DC II), Daya Bay phase I {DB Mid) and phase II {DB Full), and RENO {RN). Also 
quoted in this table are the average distances, L, to a single equivalent core (SCE) with total power 
P = X^^ii (see text for explanations, the star {*) indicates a particular treatement). The sensitivity in 
the SCE case, sin^(20i3)f;^^, is then computed to highhght possible drawback on the site configuration by 
comparison with the baseline sensitivity sin^(20i3)iini. 



a first remark, the biggest discrepancy between sin^(26'i3)iini and sin^(26'i3)fj^^ is as large as 40 % for 
the DB Mid experiment. With four to six times higher luminosity compared to DC, a farther near site 
from the cores and a closer far site, the sensitivity for the hypothetical DB Mid SCE experiment is not 
improved with respect to DC. The discrepancy between DB Mid and DB Mid SCE clearly comes from the 
wide repartition of NPP. LA cores can not be properly monitored at the DB near site. We conclude that 
the DB Mid experiment is half way between an experiment with a single far detector and an experiment 
with two identical detectors, one near and one far. One may notice that DC is an experiment where the 
reactor cores may be considered as a single equivalent core with double power since there is only a less 
than 1 % difference between DC and DC SCE. This is not the case for the RN experiment where the 
relative discrepancy between RN and RN SCE is at the level of 15 %. Daya Bay Phase II, thanks to its two 
near sites, is only marginally affected by a ~ 5 % difference between DB Full and DB Full SCE (more on 
section 18. 2p . 



8.2 The pulls analysis 

In this discussion we are interested in assessing how much a given experiment sensitivity relies on the 
knowldege of the systematics (e.^.thc number of systematics and the impact on sensitivity). The puUs- 
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approach is perfectly adapted to this analysis. The idea is to break down the total Ax^ 



Ax'(sii, 



X (-^lim) Xmin 



min X^(siim,ai,. 

{ai,...,Qjr} 



, ,a_ft:)- min x^(s, "i^ • ■ • , "a:) , 

{s,ai,...,aK] 



(where s is shorthand for sin^(2^^l3)) 



(13) 



(14) 



into sub-parts 5x^ which represent their relative contribution to the overall Ax^ of eq. ([T3|) . Since the 
sensitivity limit on sin (26*13) is computed at 90 % C.L., Ax^ has a common value for every experiment 
and is equal to 2.71. Thus, (5xf is defined as 



(i*'^ pull term) 



(15) 



with the trivial induced normalization: (5x.f = 1. 

Table fT4l shows the results of our computation for the baseline option of each experimental setup. We report, 
together with the final sensitivities discussed earlier, the relative contributions iJxf where (Jxat^ at, p are the 
observables (the contribution from the first term in eq. ([6]) for the respective detector) and the other Sxf 
are the pulls (the weight term contributions in eq. Since the role of the far detector is to determine 



Sx^ 



% 



DC I DC II DB Mid DB Full 



RN 



o 



Sx%, 



Sx' 



N2 



5x] 



3.0 % 4.3 % 



29.5 % 38.0 % 23.4 % 



1.1 % 
3.3 % 
31.2 % 



1.5 % 



34.4 % 



7i 



5X^ 
5x' 



l,abs 



scl,rel 
2 

bkg 



5x' 



Sx 
Sx 



pwr 
2 

cmp 
2 



29.1 % 

18.4 % 

6.5 % 

1.0 % 
14.7 % 
0.1 % 
0.6 % 



1.5 % 
1.3 % 
48.3 % 
1.2 % 
5.0 % 
0.8 % 
0.8 % 
0.0 % 
0.0 % 



9.0 % 

8.4 % 
6.6 % 

6.1 % 
11.8 % 
0.4 % 
27.3 % 
0.6 % 

2.2 % 



1.0 % 
0.5 % 
56.8 % 
0.1 % 
1.6 % 
0.1 % 
3.9 % 
0.1 % 
0.3 % 



7.9 % 

1.0 % 
28.5 % 
0.1 % 
0.2 % 
0.5 % 
16.9 % 

9.1 %: 

0.0 % 



'(2^13) 



13jlim 



0.054 0.028 



0.041 



0.011 



0.021 



Table 14: Relative contributions, Sxf, to the global Ax^. The higher is the Sxf contribution, the more 
the sin^(20i3) sensitivity depends on the considered parameter. In red we highlight the main contributions 
(20-60 %), in orange other significant terms (5-20 %). All the values are calculated for the base case of 
each experimental setup. For the DB Full setup, where correlations between detectors on a same site are 
possible, we take CTicI half correlated and half uncorrelated between detectors of a same site, and completely 
uncorrelated between detectors of different sites (see section [TH] for details). Note that for Daya Bay we do 
not include here batch-to-batch uncertainties. 



sin^(20i3), it is obvious that the associated residual should significantly contribute to the sensitivity. How- 
ever what we see is that all computed sensitivities mainly depend on systematics, which contribute from 60 % 
to 70 % of the overall Ax^. In the concept of identical detectors, correlated uncertainties between detectors 
should weakly impact the sensitivity (at the level of near detector "precision"). This is automatically the 
case if the near detector successfully monitors the NPP cores. Two of the quoted experimental setups reach 
this goal: DC II and DB Full. Double Chooz, with its final two detector installation, has one dominant 
systematic: the relative normalization uncertainty. The second most important contribution comes from 
the relative energy scale uncertainty. 
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Daya Bay full installation is mostly limited by the relative normalization uncertainty. The weaker impact 
of the relative energy scale uncertainties comes from the better far site distance to the NPP Cores. The 
uncertainty on the energy scale matches less the oscillation induced distortion. On the other hand, three 
of the described experimental setups still rely on theoretical knowledge of the spectrum and NPP cores 
associated uncertainties: DC I, DB Mid and RN. In particular, the DB Mid installation is sensitive to 
several systematics, especially the NPP power uncertainties. Taking data on a longer time scale (3 years, 
for instance) with such a detector configuration will not improve the sin^(20i3) sensitivity as much. Single 
core power uncertainties do not weaken the sensitivity in two particular cases: 

- if the near and far detectors have the same NPP core flux ratio contributions (this is the case of 
DC II); 

- if the far detector distance to each NPP core is the same, even if the spectra ratios are not the same in 
near/far detectors. In this particular case, the oscillation pattern is not entangled with power uncer- 
tainties in the far detector, since the Ve travelling distances are the same. Power uncertainties would 
only contribute, weakly, through the absolute normalization and correlations with other systematics 
in the near detector. 

The DB site configuration does not meet any of the above conditions. Moreover, for DB Mid, no near site 
monitors the LA I NPP. This makes the DB Mid setup an intermediate between a two identical near/far 
detector experiment and a single far detector configuration. Since the near site is farther away from the 
NPP cores (Table [T5|) . theoretical uncertainties on the spectrum have a larger impact than in DC II. Also, 
since the average distance of DB Mid far site is closer, the oscillation pattern matches slightly better with 
the energy scale associated distortion (bigger contribution than in DC II). 

The RENO far site location is the best among the quoted setups in canceling the impact of the relative 
and absolute energy scale uncertainties. However, because the site configuration does not fill any of the two 
conditions for the cancellation of the NPP core power uncertainties, this experiment relies on the precision 
with which each core power can be determined. Moreover, the near detector is a bit farther away than in 
the DC II case, which explains why the global reactor Vf, rate is less effectively determined and have a larger 
contribution than in DC II to the final sin^ (26*13) sensitivity. 

8.3 Touching the "right systematic chord" 

In the previous section we have determined the dominant systematics of each setup. In this section we 
focus on the comparison of all the reactor experiments under the assumption that systematics are known 
at the same level. Moreover, we want to illustrate the impact of the determination of the most significant 
systematics on the sensitivity of each experimental setup: 

• the single core power uncertainties; 

• the relative normalization between detectors; 

• the energy scale uncertainties (absolute and relative, between detectors). 

In the CHOOZ experiment, the power uncertainties were assessed at 0.6 % [9]. However, this estimate 
was uniquely based on the heat balance of the steam generators. Even if quite precise, this method could 
be inaccurate. Other methods, such as the external neutron flux measurements, more directly linked to 
the fission rates inside the cores, lead to a power evaluation less effectively determined with an assessed 
error around 1.5 %. This latter method allows continuous tracking of the NPP core power variations. The 
Double Chooz and Daya Bay proposals set their baseline estimates of these uncertainties at a conservative 
value of 2 %. We studied the impact of this knowledge on the sensitivity by assuming two extreme scenarios: 
in the worst case, a 3 % error, and the best case of 0.6 % precision. As a standard, we set central value 
Cpwr = 2.0 % for all the experiments. 

The relative normalization between detectors is the most significant systematic in two identical detector 
setups, with a near detector successfully monitoring the whole NPP (DC II and DB Full setups). The 
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Figure 7: Double Chooz. Daya Bay and RENO sensitivities as a function of the size of the main systematics. 
The common systematic framework is what experimentahsts believe to be achievable, without any further 
R&D. It is worth noting that the main difference between the common and the baseline cases comes from 
Double Chooz, which takes a conservative value of the relative normalization at 0.6 %. The common 
framework is used to compute the reference sin^(20i3) sensitivity of each setup (value on top of each graph). 
Then each systematic {(Jpwn o'rei, ctsci) impact on sensitivity is separately computed and illustrated as 
ratio R = sin^(26'i3)bost or worst/ sin^(26'i3)basoiinc on each graph. The overall impact changing all three 
systematics together is also illustrated with the "Total" label. Moreover we also provide a quick guess on 
sin^(20i3) sensitivity behaviour as a function of Am^^ best fit value provided by other experiments. For the 
Daya Bay Phase II experiment, where possible correlation between detectors on a same site may happen, we 
take fJroi half correlated and half uncorrelated between detectors of a same site, and completely uncorrelated 
between detectors of different sites (see section [7^ for details). 
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two available detailed quantifications of this uncertainty are the Double Chooz [51] and Daya Bay |2S] 
proposals. Each detector is estimated to measure the rate to a relative accuracy with respect to each 
other of 0.39 % (DB Full) and 0.44 % (DC II). In the Double Chooz proposal, however, a conservative 
value has been preferred for the baseline sensitivity calculation: a-^ci = 0.6 %. We will take this number 
as the worst case. The Daya Bay collaboration plans, after some R&D, to reach a relative uncertainty of 
0.2 %. This value has been taken as the best case. We set as a standard central value firoi — 0.4 % for all 
experiments. 

Even if we did not implement detector response in this simulations, we included energy scale uncertainties. 
The Double Chooz proposal quotes a^"^^ = 0.5 % and a^f^ = 0.5 %. This is taken as the common central 
values for all the experiments. However, in the CHOOZ experiment [5], the energy scale uncertainty was 
estimated at the level of 1.1 %. We thus take a 1.0 % uncertainty on both relative and absolute energy scale 
determination as the worst case. As a best case, we switched off the impact of energy scale systematics on 
the sensitivity. 

Although the impact of the Am§]^ uncertainty on the sensitivity to sin^(20i3) is negligible with the current 
knowledge, the central best fit value on this parameter from other experiments such as MINOS [H], K2K [7] 
and SuperK may have some influence on the reactor experiment sensitivities. Wc thus include the current 
Am§]^ bounds on the sensitivity computations: the "worst case" is taken to be Ato|]^ = 2.0 lO"'^ eV^ and the 
"best" one Am|j^ = 3.0 10^'^ eV^ (2 a bounds from [T^). The standard central value for all the experiments 
has been taken to be Amgj = 2.5 10~^ eV^. 

We illustrate in Figure [7] these three systematic scenarii: best, central and worst cases. Each contribution 
is assessed separately, but we also show on this graph the total impact by summing the effect of the 
three discussed systematics (cTpwr, Crci and CTsci) to their respective best, central and worst values. We also 
illustrate the impact of the true central value of Am|]^ on the sensitivity. In this representation, we show the 
ratio of the computed sensitivity sin^(20i3)b,w for the best (resp. worst) case over the sensitivity sin^(20i3)c 
for standard central systematic values. The "Total" bar shows that in Daya Bay and RENO the sensitivity 
can vary from 0.6 to 1.2-1.3 of the baseline case, for experimental systematics ranging from a "best" to 
a "worst" scenario. In the case of Double Chooz, the impact of systematics is less significant, at the level 
of 20 % on both sides. In the Double Chooz and DB Mid cases, the sensitivity could be worsened for best 
fit values of Ato§j below 2.5 10"'^ eV^. In DB Full, the sensitivity is quite stable over the current allowed 
range for Ato|j with only a 5 % effect on sensitivity. 

Double Chooz is an optimized experiment in the sense of robustness with respect to systematics for a goal 
sensitivity in the 0.02-0.03 range. Daya Bay phase II is adequate to reach a sensitivity at the level of 0.01. 
However, a simpler experiment for this class of sensitivity would be a scaled-up variant of Double Chooz, 
with a very close near site at 150 m and a 1.5 km baseline for the far site. At the Diablo Canyon power 
plant [Hj, where two 3.19 GWth cores are operational and modest civil engineering works would be required, 
four 20 t detectors (2 near, 2 far) would give a sensitivity of 0.013 after 3 years of data taking. 

9 Conclusion 

In this work we have presented a detailed comparative analysis of the sensitivities to sin^ (26*13) of upcoming 
and proposed reactor experiments. 

We have first calculated the sensitivities using all available data published by the respective collaborations 
for the baseline of both the systematical uncertainties and the experimental setups. Our results are generally 
in good agreement with the sensitivities quoted by the collaborations: 0.054 for Double Chooz Phase I; 0.028 
for Double Chooz Phase II; 0.041 for Daya Bay "mid"; 0.0089 for Daya Bay "full"; and 0.021 for RENO. 
In the case of Daya Bay, we have additionally evaluated the impact of the proposed filling procedure, 
with pairs of near-far detectors filled with the same scintillator batch and 4 different scintillator batches. 
If the hydrogen mass fraction is controlled to 0.5 % between different batches, the sensitivity worsens 
slightly to 0.0093. We also examined, still for Day Bay, the until-now implicit assumption that the errors 
of the relative normalizations between the 8 detectors are fully uncorrclatcd. In reality this is the most 
optimistic scenario, since a part of these uncertainties may come from site-dependent systematics. In the 
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most pessimistic scenario (all relative normalizations fully correlated for the same site), the Daya Bay 
sensitivity would be 0.012. 

An important result of this work is that the total thermal power available for an experiment, a figure of 
merit that has been often used as a strong argument to "rank" different projects, has a modest impact on 
the success of an experiment. Large powers are only available in multi-core sites, which are very difficult to 
monitor. The associated systematics can be overwhelming with respect to the benefit from the statistics. 
This is very nicely exemplified by the case of RENO, which would reach the same sensitivity with just 2 of 
its 6 reactors on, and by Daya Bay "mid" , which results to be just half way between Double Chooz Phase I 
and Double Chooz Phase II. Wc have illustrated how optimal is the use of the available thermal power in 
each site through a comparison of the real experiments with ideal "single core equivalent" setups: Double 
Chooz and Daya Bay "full" arc nearly optimal, Daya Bay "mid" and RENO do not make full advantage of 
their huge power. Daya Bay pays for the complexity of its nuclear power plant by the inevitable construction 
of two near sites. 

Wc have carried out a detailed and unified x^— pull analysis for all the experiments, under common assump- 
tions for all systematic uncertainties. This study allowed us to compare all projects on an equal footing and 
evaluating the impact of each single systematic uncertainty on the sensitivity of each experiment. In this 
calculation the Double Chooz baseline sensitivity is 0.0235 and a single systematic dominates, which is the 
error of the relative normalization between near and far detector. This shows that taking into consideration 
its small mass. Double Chooz is an optimized experiment. With respect to the other projects, the Double 
Chooz sensitivity has a more pronounced dependence on the best fit Am^g, with sin^(20i3)iini ranging from 
0.020 to 0.031 for the currently 2 cr— allowed Arnf^ interval. 

Daya Bay "full" proves to be robust with respect to the size of the systematical errors and to variations of 
Am^3 and is the only partly approved project with an achievable sensitivity potential below 0.01. The Daya 
Bay sensitivity is dominated by the accuracy of the relative normalization between detectors and the degree 
of correlation existing on a same site between the latter parameters. With the knowledge of the systematics 
we have today, the Daya Bay sensitivity would be on average between 0.0094 and 0.0123, depending on the 
degree of correlation of the systematics between detectors on the same site. 

Contrary to Double Chooz and Daya Bay, the sensitivity of RENO is largely degraded by the uncertainties 
of the reactor powers and fuel composition. This, again, shows that the site is not optimal. Nevertheless, 
due to the large target mass and optimal baseline, a very competitive sensitivity, around 0.02, is achievable. 

This pull analysis also shows that the impact of the backgrounds on the is minor with respect to other 
systematics. Backgrounds, at least in our simulation, are therefore not critical in any of the analyzed 
experiments. 

Taken into consideration all the above results, we come to a conclusion about the features of the optimal 
experiment to approach a sin^(20i3) sensitivity of 0.01. Where by "optimal" we mean a robust sensitivity, 
the simplest configuration, the minimal amount of civil works and the smallest mass. Such an optimal 
experiment would have a nuclear power plant layout as simple and powerful as Double Chooz; a favorable 
topography for sufficiently deep underground laboratories, a very close, single near site, a ~ 1.5 km baseline 
for the far site and a total target mass of about the half of Daya Bay. Diablo Canyon, California, is a good 
example of an suitable site for such an experiment. 
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Appendix 



In this section we first describe the way we computed event rates, and then the analysis implementation 
with detailed systematic inclusions. 



Event rates 

The visible energy inside the detector has a simple expression as a function of the positron energy or neutrino 
energy of the inverse /3-decay reaction: 

Svis = E,+ + ~E,~ E'J"' + 2 me , (16) 

where E^'^ = Af„ — Mp + rrie is the Ve threshold energy of the reaction. The event rates produced in reactor 
R and recorded in detector D per visible energy bin [Ef. -B^+i] may be written as 

iVf'^(0i3,Am2i)= / AL AE,^ AE, —^{E,+ , E,) (17) 



with 



'R,D 



——(£;,+ ,£,) = ^l^h{LR,D,L)a{E,)^^^''{E,)e{E,+ ) (18) 



xR{E,+ ,E,)P,e{E,,L,9i3,^ml^) , 

where h^Lji Oi L)^ a{E^), e{E^+), 4'^'^{E^), R{E^+ ,E^) stand for the finite size effect function the 77e 
inverse /3-decay reaction cross section, the detection efficiency, the 77e fiux from reactor R in detector D 
and the energy response of the detector respectively. The generic normalization factor M^'^ is the product 
of the experiment life time by the number of available free protons inside the target volume, the global load 
factor of reactor R and the dead time of detector D. The flux from reactor R in detector D is described 
in term of the isotope composition as: 

r<R 

^«'^(i?.) = P^E^'^f''(^^) (19) 

where £ — ^^^U, ^^^U, ^"^^Pu, ^"'^Pu labels the most important isotopes contributing to the flux, Cf^ 
is the relative contribution of the isotope £ to the total reactor power [C^ = Nf^Ef^/P^, and A^^** is the 
number of fissions per second of isotope i), Ef^ is the energy release per fission for isotope £ and is the 
thermal power of reactor R. In eq. (|19p . (l)f'^{E^) is the energy differential number of neutrinos emitted 
per fission by the isotope £, and we adopt the parameterization for the (j)f'^{E^) from ref. jTB]. For the Cf^ 
we take a typical isotope composition in a nuclear reactor given in eq. P|). 

We take for Pee the oscillation probability expressed in eq. Wc assume in this article a constant efficiency 
of e(-E'e+) = 80 % and a Gaussian energy response function 




with an energy resolution of p{Eg+) = 8 %/ \/E^+[MeV]. The finite size effect function h{LR^£,,L) is 
assumed to be a Dirac distribution (pointlike sources and detectors) in this paper. The total event rates in 
i**^ bin of detector D is then simply expressed as 

iVf(0i3,Am2,)= ^ iV«'^(0i3,Amii). (21) 

R=Ri....Mn, 
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analysis and systematics inclusion 



The computed event rates, N^'^ , are then mcluded m a pull-approach analysis 
between the systematic uncertainties are properly included: 

/ K 



X 



mm 



J2 + J2 a,W,-l, 



aw 



i=l,. 

D=Di, 



..Dr. 



fe = l 



fe,fc' 



where correlations 



(22) 



with, 



and Uf are given by eq. 



The Q!fe and SR. coefficients are described in Table fT51 



(23) 



The SPj, coefficient represents the shift of the i*"^ bin of detector D spectrum due to a 1 cr variation in the 



k*'^ systematic uncertainty parameter ak ■ Most of the systematics are expressed as function of ' of Bf 
quantities already described previously. The energy scale systematic coefficients in Sf^. are defined through 
MP which follows the relation 
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R=Ri,...,Ri. 
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R.D 



where 



/ dL / d^,+ / 

^0 J Ei-m^ Je 



r. + oo d'^M^'^ 



-{Ee+,E^) . 



(24) 



(25) 



It is often assumed in the pull-approach that Wk,k' = <5fc,fc'- If we keep this definition of Wk,k' then we 
are faced with the problem that the reactor spectrum shape uncertainties may contribute to the absolute 
normalization error and the fuel composition uncertainties may contribute to the reactor core power errors. 
If we want to get rid of these free contributions we could use two methods: 

1. the first one consists to infer that Qfcmp/C/^ = 

2. the second one is to introduce additional weight terms in the definition (|22p : Oicmp.eCf /ecmp 

The first method allows disentangling fuel composition from power uncertainties. However, this assumption 
is a bit too restrictive since in practice the sum ctcm^/Cf^P^ = is only constrained at the knowledge 
level of the power of reactor R. The second method has the advantage of allowing an estimate of the level 
of contribution of this systematic to the power uncertainties. Moreover this simply leads to redefining the 
Wk,k' matrix as 

r<RrR 



with 



ej',k,k' 2 
.R. ^cmp 



if fc, k' e {fc^ + 1 



jcmp.R 



k^' = 5A^d 







"0 ~ 
otherwise 



, fc(f + 5}, 
+ iVb + A^r + 2 + 4{R - 1) 



(26) 



(27) 



With these definitions, ffcmp determines at which level the fuel composition uncertainties are allowed to 
contribute to the core power errors. One wants typically that fuel composition uncertainties contribute 
within the allowed region of power uncertainty. Thus, 



-cmp 



^ pwr 



Regarding the fuel composition uncertainties, CTcmp may be assessed roughly at the level of cTpwr 

aLr. = 2-3%. 



(28) 

uncertainties: 
(29) 
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Relative Energy scale {Ng = 


Nd + 2) 








in Di 


A^.H 


- 1 




"scl 


in Dn^ 


A^s + 


Ad 






Backgrounds (N^ = 2iVd + 2) 








accidentals in Di 


A,H 


- 1 


Sl 1,4 




accidentals in D^^ 


A^s + 


Ad 






cofainogemco m ui 


A, + A^d + 1 


B2 2,4 


"-B2 


cosmogcnics in Dn^ 


A^s + 


2Ard 






proton recoils in Di 


Ns + 2Ad + 1 






proton recoils in -Dat^ 


Ns + 


SA^d 




B3 


Reactor spectrum shape {Ng 


= 5Ad- 


+ 2) 






in bin 1 


A,H 


- 1 


CTshpAf 


ashp.i 


in bin A'^b 


A^s + 


Ab 




Q!shp,Arb 


Reactor power (Ag = 5A^d + 


A^b + 2) 








from i?i 


A^sH 


- 1 


'-^ pwr * i 


^pwr 


from Rn^ 


Afs + 


A^r 


(Jpwr 




Reactor R composition (A"s = 


= 5A^d + 




A^r + 2 + 4(i? - 


1)) 


from 235u 


AfsH 


- 1 


CcmpA^,- 




from Pu 


A^sH 


-2 


fcmpjV,; 


^cmp,i? 


from 238 U 


A^sH 


-3 


CcmpA^i 


238y 


from ^''^Pu 


A^sH 


-4 


Ccmp-'V,; 





Table 15: Systematic parameters table. We used the following definitions: Ab is the number of bins in the 
reactor spectra, Ad is the number of detectors in the experiment, N,- is the number of reactor cores in the 
NPP, A^s is short for the number of previously defined systematic parameters. For specific values of Cabsj 
Crei, CTghp, o'qu, (Jpwr, cTcmp, we refer to the experiment comparison sections in]-[Sl 



27 



References 



B. T. Cleveland et at, Astrophys. J. 496 (1998) 505; J.N. Abdurashitov et al. [SAGE Collaboration], 
J. Exp. Theor. Phys. 95 (2002) 181, |aitro -ph/0204245; T. Kirsten et al. [GALLEX and GNO Collab- 
orations], Nucl. Phys. B (Proc. Suppl.) 118 (2003) 33; C. Cattadori, Talk given at Neutrino04, June 
14-19, 2004, Paris, France. 



Super-K Collaboration, S. Fukuda et al, Phys. Lett. B 539 (2002) 179, hep-ex/0205075 J. Hosaka et 



al. |hep-ex/0508053 



SNO Collaboration, Q.R. Ahmad et al, Phys. Rev. Lett. 89, 01 1302 (2002), |nucl-ex/0204009[ 



B. Aharmim et al., Phys. Rev. C 72, 055502 (2005), nucl-ex/0502021 



Super-K Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 81 (1998) 1562, hep-ex/9807003 Y. Ashie 



et al, Phys. Rev. D 71 (2005) 112005, hep-ex/0501064 



K. Eguchi et al [KamLAND Collaboration], Phys. Rev. Lett. 90 (2003) 021802, hep-ex/0212021 



T. Araki et al [KamLAND Collaboration], Phys. Rev. Lett. 94, 081801 (2005), hep-ex/0406035 



K2K Collaboration, E. Aliu et al, Phys. Rev. Lett. 94, 081802 (2005) , |hep-ex/041 1038| M. H. Ahn et 
aZ.,|hep-ex/0606032, 



B. J. Rebel [MINOS Collaboration], arXiv:hep-ex/0701049 



M. ApoUonio et al. [CHOOZ Collaboration], Phys. Lett. B 466, 415 (1999), [hep"-ex/9907037( Eur. 
Phys. J. C 27, 331 (2003), .hep-ex/0301017 , 

G. Ahmonti et al. [BOREXINO Collaboration], Astropart. Phys. 8, 141 (1998). 

B. Pontecorvo, J. Exptl. Theoret. Phys. 33 (1957) 549 [Sov. Phys. JETP 6 (1958) 429]; J. Exptl. 
Theoret. Phys. 34 (1958) 247 [Sov. Phys. JETP 7 (1958) 172]; Z. Maki, M. Nakagawa and S. Sakata, 
Prog. Theor. Phys. 28 (1962) 870. 



M. Maltoni et al. New J. Phys. 6, 122 (2004), hep-ph/0405172 



T. Schwetz, Phys. Scripta T127 (2006) 1, | hep-p h/0606060[ 

G. L. Fogh, E. Lisi, A. Marrone, D. Montanino and A. Palazzo, Phys. Rev. D 66, 053010 (2002), 



hep-ph/0206162 



Schreckenbach K. et al, Phys. Lett. B160,325(1985); Hahn A.A. et al, Phys. Lett. B218, 365(1989); 
F. Von Feilitzsch et al Phys. Lett. B118 , 162(1982). 

P. Vogel and J. Engel, Phys. Rev. D 39 (1989) 3378. 



C. Bemporad, G. Gratta and P. Vogel, Rev. Mod. Phys. 74 (2002) 297, hep-ph/0107277 



P. Huber and T. Schwetz, Phys. Rev. D 70, 053011 (2004) |hep-ph/0407026 
K. Anderson et a?., | hep-ex/0402041 



P.Antonioh et al, Astropart. Phys.7, 4 (1997), 357-368 

B.P. Heisinger, Muon-induzierte Produktion von Radionukliden, dissertation at the Physik Department 
E15, TU-Miinchen, (1998). 

A.Tang, G. Horton-Smith, V. A. Kudryavtsev, A. Tonazzo, Journal-ref: Phys.Rev. D74 (2006) 053007. 

T. Hagner, R. von Hentig, B. Heisinger, L. Oberauer, S. Schonert, F. von Feilitzsch, and E. Nolte, 
Astropart. Phys. 14, 33 (2000). 



28 



[24] 
[25] 

[26] 
[27] 
[28] 
[29] 
[30] 



F. Ardcllicr et ai, h ep-ex/0405032[ S. Berridge et ai, |hep-ex/0410081[ M. Goodman, Th. Lasserre et 
g/., hep-ex/0606025. 



X. 



Guo 



hep-ex/0610024 



_[Daya Bay Collaboration], "hep-ex/0701029' Yifang Wang, 

Seminar of D. Jaffc at CEA/Saclay on October 23 2006, 



www-dapnia. cea.fr/Phocea/file .php?class=std&f ile=Seminaires/ 1425/ 1425 .pdf 
http : //neutrino . snu . ac .kr/RENO/ 

http : //www .ba. inf n. it/~now/now2006/Progrciin/program.html 



M. Aoki et aL, |hep-ex/0607013 
J.C. Anjos et 



hep-ex/0511059 



CEA, Nuclear power plants in the world (2005). 

http : //132 . 166 . 172 . 2/f r/Publications/tr ilogie/Elecnuc-2005 . pdf 

CEA, Energy handbook (2005). 

http : //132 . 166 . 172 . 2/f r/Publications/trilogie/Inf os-Energie-2005 . pdf 



Achkar B. et ai, Phys. Lett. B374, 243(1996). 
H. Minakata, et ai, Phys.Rcv. D68 (2003) 033017, |hep-ph/02111H" 



[31] 
[32] 

[33] P. Huber, et ai, Nucl.Phys. B665 (2003) 487-519, |hep-ph/0303232 
[34] G. Mention, PhD. Thesis, [http:/ /doublechooz.in2p3.fr7 

[35] 
[36] 



M. ApoUonio, et al, Eur. Phys. J. C27, (2003) 331, |hep-ex/0301017 
S. Eidelman et ai, Phys. Lett. B592 (2004). 



29 



